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ABSTRACT 

We present theoretical models for the evolution of T Tauri stars surrounded by cir- 
cumstellar discs. The models include the effects of prc-main-sequence stellar and time 
dependent disc evolution, and incorporate the effects of stellar magnetic fields acting 
on the inner disc. For single stars, consistency with observations in Taurus-Auriga 
demands that disc dispersal occurs rapidly, on much less than the viscous timescale 
of the disc, at roughly the epoch when heating by stellar radiation first dominates 
over internal viscous dissipation. Applying the models to close binaries, we find that 
because the initial conditions for discs in binaries are uncertain, studies of extreme 
mass ratio systems are required to provide a stringent test of theoretical disc evolution 
models. We also note that no correlation of the infra-red colours of T Tauri stars with 
their rotation rate is observed, in apparent contradiction to the predictions of simple 
magnetospheric accretion models. 

Key words: stars: pre-main-scquence stars: magnetic fields - accretion, accretion 
discs - binaries: general - stars: rotation - circumstellar matter 



1 INTRODUCTION 



Observations of pre-main-sequence T Tauri stars suggest 
that stellar magnetic fields are often important in the inner 
regions of surrounding circumstellar discs. The field acts to 
provide a link between the star and its disc, and the result- 
ing torques lead to truncation of the disc interior to some 
magnetospheric radius 7? m (see e.g. Pringle & Rees 1972; 
Ghosh & Lamb 1979). This divides the accretion flow into 
two regimes; an outer disc where inflow is slow and the sur- 
face density of the disc material relatively large, and an inner 
evacuated 'hole' where the magnetic field is dominant and 
infall onto the star proceeds at approximately the free-fall 
velocity. 

Theoretical work shows that stellar magnetic fields of 
the order of ~ 1 kG are a prerequisite for this mode of ac- 
cretion (Konigl 1991). Magnetic fields of this strength have 
been directly detected in a few weak-line T Tauri stars 
(WTTS - Basri, Marcy & Valenti 1992; Guenther & Emer- 
son 1995, 1996), while the high X-ray flux (Montmerle 1992; 
Neuhauser et al. 1995) and flaring activity (Preibisch, Zin- 
necker & Schmitt 1993) seen in the general T Tauri popula- 
tion suggest that strong magnetic fields are common. Several 
consequences of magnetospheric accretion are also observed, 
including mass infall at free fall velocities (Calvet & Hart- 
mann 1992; Edwards et al. 1994; Martin 1997; Johns- Krull & 
Hatzes 1997), the presence of hotspots in classical T Tauri 
stars (CTTS - Bouvier et al. 1993), and infra-red colours 



that are too red to match disc models that extend to the 
stellar surface (Bertout, Basri & Bouvier 1988). 

Theoretical models that aim to describe the star-disc in- 
teraction are abundant (e.g. Ghosh & Lamb 1979; Campbell 
1992; Yi 1994; Lynden-Bell & Boily 1994; Wang 1995, 1996; 
Lovelace, Romanova & Bisnovatyi-Kogan 1995; Ostriker & 
Shu 1995). There are important differences between these 
models, but a generic prediction is that the magnetospheric 
radius should lie close to the corotation radius (R c , where 
the the disc material is stationary in the rotating stellar 
frame), if the magnetic linkage is to be effective in regulat- 
ing the star's angular momentum (Wang 1995; Armitage & 
Clarke 1996). 

In recent studies, Kenyon, Yi & Hartmann (1996) and 
Meyer, Calvet & Hillenbrand (1997), compared the infra-red 
colours of a sample of classical T Tauri stars with the predic- 
tions of several steady-state magnetic and non-magnetic disc 
models. They concluded that the magnetically disrupted 
discs provided a better fit to the observed colours, and were 
able to place some constraints on the location of the mag- 
netosphere relative to corotation. 

In this paper, we employ time-dependent theoretical 
models to investigate how magnetically truncated discs 
evolve with time. Because the viscous timescale in the region 
where infrared emission originates (R < 1 a.u.) is small, we 
expect that steady-state models such as those used previ- 
ously will suffice for calculating the infrared colours (assum- 
ing that the inner regions are stable). However the use of a 
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Figure 1. HR diagram for the stellar models used in this paper, 
with stellar masses Af* = 0.9885, 0.4375, 0.1901, 0.10 and 0.057 
Mq. Tracks are plotted started (arbitrarily) where = 3 -Rq, 
with isochrones shown relative to this time at At = 1, 2, 4, 8 Myr. 

time- dependent approach permits us to study how rapidly 
the system evolves in, for example, the infrared colour-colour 
plane. 

Section 2 describes the numerical models used in this 
work, which are applied in Section 3 to follow the evolution 
in the infra-red colour-colour diagram of discs surrounding 
magnetic T Tauri stars. For a limited number of systems 
additional information, in the form of measured stellar ro- 
tation periods, is available, and we consider the constraints 
this imposes in Section 4. Section 5 applies our models to 
discs in close binary systems. Section 6 summarises our con- 
clusions. 



2 NUMERICAL MODELLING 

The numerical models employed in this work are derived 
from those used previously by Armitage & Clarke (1996) 
for modelling T Tauri systems with discs. We follow the 
evolution of a theoretical pre-main-sequence stellar model 
surrounded by a geometrically thin a-prescription (Shakura 
& Sunyaev 1973) accretion disc. A treatment that fully cou- 
pled the stellar evolution calculation to the disc would be 
essential for the early stages of accretion, here we only con- 
sider accretion rates low enough that the stellar evolution 
and accretion problems can be considered separately. 

2.1 Stellar model 

Figure [j] shows the Hertzprung-Russell diagram for the the- 
oretical pre-main-sequence models employed in this paper. 
We construct our stellar models using the most recent ver- 
sion of the Eggleton evolution program (Eggleton 1971, 
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Figure 2. Time evolution of the stellar radii, for models with 
stellar masses (from top downwards) M* = 0.9885, 0.4375, 0.1901, 
0.10 and 0.057 Mq. 



1972, 1973). The equation of state, which includes molecu- 
lar hydrogen, pressure ionization and coulomb interactions, 
is discussed by Pols et al. (1995). The initial composi- 
tion is supposed to be uniform with a hydrogen abundance 
X = 0.7, helium Y = 0.28, deuterium A D = 3.5 x 10" 5 
and metals Z — 0.02 with the meteoritic mixture deter- 
mined by Anders and Grevesse (1989). Hydrogen burning 
is allowed by the pp chain and the CNO cycles. Deuterium 
burning is explicitly included at temperatures too low for 
the pp chain. Once the pp chain is active hydrogen is as- 
sumed to burn to He 4 via deuterium and He 3 in equilibrium. 
The burning of He 3 is not explicitly followed. Opacity tables 
are those calculated by Iglesias, Rogers and Wilson (1992) 
and Alexander and Ferguson (1994). An Eddington approx- 
imation (Woolley and Stibbs 1953) is used for the surface 
boundary conditions at an optical depth of r = 2/3 so that 
low-temperature atmospheres, in which convection extends 
out as far as r « 0.01 (Baraffe et al. 1995), are not modelled 
perfectly. However the effect on observable quantities (see 
Kroupa and Tout 1997) is not significant in this work. 

The tracks in figure 1 are for four pre-main-sequence 
stars of masses 0.9885, 0.4375, 0.1901, 0.1 and 0.057 M© 
shown from a radius of 3 Rq until they begin hydrogen 
burning at the main sequence. The most massive model 
establishes a radiative core well above the main sequence 
and so follows a Henyey track in its last stages of contrac- 
tion. The other models remain fully convective although the 
0.4375 Mq star will develop a radiative core on the main se- 
quence. The actual starting point for all of our models was 
an accreting protostellar track in which material of the zero- 
age composition and with the surface state is accreted at a 
rate of 10 M©yr _1 . Accretion was interrupted when the 
various total masses were reached and the stars allowed to 
evolve freely with normal photospheric boundary conditions 
and constant total mass. All these models have relaxed from 
their accreting state by the time they reach 3 Rq and are 
by that time descending Hayashi tracks as typical pre-main- 
sequence stars. Figure 2 illustrates the evolution of radius 
with time for each of these stars. The three lower-mass ob- 
jects show a plateau lasting between one and five million 
years. This corresponds to the delay in shrinkage while the 
stars burn their primordial deuterium. The duration of this 
plateau will depend on the actual value assumed for the pri- 
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mordial deuterium abundance. The two more massive stars 
have already burnt all their deuterium by the time they have 
shrunk to 3 Rq because their cores are correspondingly hot- 
ter. 

The early accretion history of protostars has been 
considered by Stahler (1988), and by Hartmann, Kenyon 
& Cassen (1997). How the accretion influences pre- main- 
sequence evolution depends upon whether accretion at early 
times is via a magnetospheric or boundary layer mode, and 
on the amount of mass gained during high accretion rate out- 
bursts. Here we are primarily interested in the later phases 
of accretion, for which purpose a rather arbitrary assignment 
of t — suffices. We choose the time at which _R* = 3 Rq 
because it corresponds to roughly to the upper limit on the 
inferred stellar radii of CTTS (Hartigan, Edwards & Ghan- 
dour 1995), and the radii used for models of stars suffering 
FU Orionis outbursts at higher accretion rates than those 
modelled here (Bell & Lin 1994). 



2.2 Disc model 

We consider a geometrically thin non-self-gravitating disc 
evolving under the influence of viscous and stellar magnetic 
torques. The evolution of the disc surface density distribu- 
tion ~E(R,t) is described by (e.g. Livio & Pringle 1992), 



as 
at 
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Here v is the kinematic viscosity, M the stellar mass, and 
B z and B$ the vertical and azimuthal components of the 
stellar magnetic field at the surface of the disc. We make the 
standard assumption that the vertical field B z has a dipolar 
fall-off with radius at the magnetospheric radius (generally 
at several R* ) , and use the expression for the azimuthal field 
given by Wang (1995), 
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where Q(R) is the Keplerian velocity at radius R in the 
disc and fi, the angular velocity of the stellar surface. In 
protostellar discs where the Shakura-Sunyaev a parameter 
is probably <JC 1, this expression is likely to be valid if rapid 
reconnection in the corona limits the growth of B$. Different 
field geometries arise if \B^,/B Z \ is able to grow to a large 
value (Bardou & Heyvaerts 1996). 

The kinematic viscosity is described using solutions to 
the vertically averaged equations for thin discs given in 
Faulkner, Lin & Papaloizou (1983). The analytic opacity 
fits used by Bell & Lin (1994) are used, so that in any 
one opacity regime the solution for i/(a, R, S) is a power- 
law in surface density, radius, and a parameter. Transitions 
to different dominant opacities occur as the disc becomes 
cooler at greater radii and later times. At an accretion rate 
of M = 10~ 8 M Q yr -1 , and radii between 5 and 80 Rq, 
these solutions agree with the detailed vertical structure in- 
tegrations presented by Bell & Lin (1994), using the same 
opacities, to within typically 20-30 %. 

At low accretion rates, the thermal structure of the disc 
will be determined primarily by the heating from the central 



star, rather than viscous dissipation within the disc. When 
the accretion rate falls low enough that this occurs, we keep 
v(R) constant for the remainder of the evolution. 



2.3 Spectral energy distributions 

We calculate the disc spectral energy distribution assum- 
ing that each annulus in the disc radiates locally the energy 
input from viscous heating and reprocessing of stellar ra- 
diation. For models including a stellar magnetic field the 
work done on the disc material by magnetic torques is also 
included. The total heating of one side of the disc is then, 

Q + = Qt + Qt + Qt, (3) 
where the viscous term is the usual expression, 

Qt = |i/En 2 , (4) 
and the magnetic contribution is given by 

Qb = i(£*B«)*=Hii 2 |n-a.|. (5) 

The component Qp~ due to reprocessing of stellar radiation is 
computed using the standard expressions given by Kenyon & 
Hartmann (1987; see also Adams & Shu 1986). We assume 
a flat disc, a stellar spectrum given by a Planck function 
B v (T e ), and allow for the partial absorption when the disc 
starts to become optically thin. 

Treating the disc as an isothermal slab, the energy loss 
rate per unit area through one side of the disc Q~ is, 



Q~ =nB v (T disc )[l-2E 3 {T )] 
where, 



£3(1-0) 
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(7) 



To(y) = k(^)E is the optical depth normal to the disk plane, 
and Tdisc is determined by the requirement that Q + — Q~ . 

The dust opacity k{v) is that used by Menshchikov & 
Henning (1997) in their modelling of the source L1551 IRS 
5, and incorporates a mixture of carbon and silicate com- 
ponents. Dust opacities vary considerably with the adopted 
grain model (Menshchikov, private communication, see also 
figure 15 of Menshchikov & Henning 1997), and this will 
modify the SED as the disc becomes optically thin. Addi- 
tionally in the inner regions of the disc, where T > 1500 K, 
the dust is destroyed and the remaining opacity, owing to 
molecules, is greatly reduced (Calvet, Patino & Magris 1991; 
Alexander, Johnston & Rypma 1983). This is not included 
in our current calculations. However, at low accretion rates, 
when the temperature distribution in the disc is dominated 
by the reprocessing of stellar radiation, we expect this to 
be an important complication only in the innermost regions 
of the disc. For example, for a star of effective temperature 
T» = 4000 K surrounded by a flat disc, the temperature falls 
below 1500 K within 3-R*. In our evolutionary models the 
magnetosphere generally extends to at least this radius at 
late times, and so we do not expect large changes in our 
results if dust destruction was treated consistently. 

For each model, we compute the spectral energy dis- 
tribution (SED) for two system inclinations, a face-on disc 
system and one inclined at 60°. Infra-red colours are then 



© 0000 RAS, MNRAS 000, 000-000 



4 P.J. Armitage, C.J. Clarke and C.A. Tout 



computed using the observed colours of main- sequence stars 
of given effective temperature (Kenyon & Hartmann, 1995), 
to which are added the contribution from the disc using the 
zero points in the UKIRT system (e.g. as quoted in the HST 
NICMOS reference manual). 



3 COLOUR EVOLUTION OF 

MAGNETOSPHERIC ACCRETION MODELS 

3.1 Initial conditions 

The basic properties of T Tauri discs are now known obser- 
vationally. Typical disc masses are ~ 1(T 2 M Q (Osterloh & 
Beckwith 1995), typical accretion rates are ~ 10 -8 M©yr _1 
(Gullbring et al. 1998), and typical disc lifetimes are a 
few Myr (Strom 1995). We aim to set up initial conditions 
that are simultaneously compatible with these observations, 
while noting that observationally there is considerable scat- 
ter in all of these quantities. We especially caution as to the 
possible systematic errors in estimating ages from pre-main- 
sequence tracks (Tout, Livio & Bonnell 1998). 

For our magnetic models we use the M — 0.4375 Mq 
stellar model, and an initial steady-state disc with an ac- 
cretion rate of 3 x 10~ 7 M©yr _1 , at the upper end of the 
T Tauri accretion rate distribution. The initial disc mass is 
taken as 0.1 M*. With these parameters, experimentation 
showed that an a of 0.02 produced sensible values for the 
disc lifetime (in the Myr range) and initial radius (~ 30 a.u., 
though the disc would expand as it evolves). This value of a 
is also consistent with the findings of Hartmann et al. (1998), 
who considered the evolution of rather different self-similar 
disc models. These runs have a zero-torque outer bound- 
ary condition imposed at 30 a.u., and use 800 radial mesh 
points. 

The models also require specification of the stellar mag- 
netic field and stellar rotation period. We parameterise the 
stellar field as, 



B* = Bi 



(8) 



and take Bo between 250 G and 2 kG. For the stellar rotation 
rate, we assume that P, varies smoothly between an initial 
spin period Pi = 7d, and a final spin period Pf — 3.5 d, con- 
sistent with observations (Bouvier et al. 1995). The transi- 
tion is smoothed according to, 

(Pi ~ Pi) 



P*=Pt + 



1 + (t/twTTs) 2 



(9) 



where the time taken to reach a WTTS spin period, iwTTS, 
is taken to be 3 Myr. 

3.2 Results 

Figures ^ and ^ shows the model tracks in the (K-L)-(K- 
N) colour plane, together with evolution of the colours with 
time. Also plotted is the observed distribution of Class II and 
Class III sources in Taurus- Auriga. The data are taken from 
the compilation of Kenyon & Hartmann (1995) and are not 
dereddened. The stars comprising the sample have a fairly 
broad distribution of ages (Kenyon & Hartmann 1995; see 
especially Fig. 16), with a typical age of ~ 1 Myr. A 'gap' in 
K-N between 1 < (K — N) < 2 is clearly seen, for this sample 
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Figure 3. Tracks in the K-L, K-N plane for magnetically trun- 
cated disc models surrounding an M* = 0.4375 Mq star, with 
B = 250, 500, 1000, and 2000 G respectively (higher B models 
are redder in K-N). The dashed lines are calculated for a face-on 
disc system, the solid lines for a disk inclined at 60°. Symbols 
are from the compilation of observations by Kenyon & Hartmann 
(1995). 
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Figure 4. Time evolution of the infra-red colours corresponding 
to the tracks shown in Fig. 3. The dashed lines are calculated for 
a face-on disc system, the solid lines for a disk inclined at 60°. 
Note that the absolute time scaling is entirely determined by the 
assumed value of a, which is here chosen as a = 0.02 in order to 
roughly match the observed lifetime of T Tauri discs. 
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the only systems lying in the gap are three binaries without 
separate colour information for the two components. 

For a given K-L, models with stronger stellar magnetic 
fields are generally redder in K-N. Those models with lower 
Bo , where for high accretion rates i? m is significantly smaller 
than R c , seem to fit the envelope of the CTTS observations 
better than models with higher stellar fields, as noted pre- 
viously (Kenyon, Yi & Hartmann 1996). 

Our calculation of disc colours is crude - discs are prob- 
ably not flat, and their spectra in the infra-red is likely to de- 
viate substantially from a blackbody. Thus the tracks taken 
by the models are more significant than the absolute match 
with the data. Two phases of evolution occur. In the first 
phase, the K-L colour declines smoothly to photospheric val- 
ues, while the K-N colour remains red. This first phase is 
driven by the decline in the accretion rate, which increases 
the relative importance of magnetic torques compared to the 
viscous torques in the inner disc, and moves the magneto- 
sphere outward towards corotation, destroying the L flux 
from the disc and reducing the K-L colour. The N flux, 
sampling cooler parts of the disc unaffected by the mag- 
netosphere, remains strong. This longer wavelength emis- 
sion lasts until the inner disc becomes optically thin. The 
inclusion of the magnetic torques does not assist in pro- 
ducing a sharp transition between disc-like and stellar K-N 
colours, which instead proceeds on the slow viscous timescale 
of the outer disc (note that Fig. ^ is plotted on a logarithmic 
timescale, and thus the models predict that the amount of 
time spent in the unpopulated gap region is roughly equal to 
that spent in the allowed CTTS part of the diagram). The 
observation of a 'gap' thus unequivocally requires disc clear- 
ing to take place on a timescale much shorter than expected 
from purely viscous disc processes, even if stellar magnetic 
fields are strong enough to rapidly evacuate the inner re- 
gions of the disc out to corotation. Since mid-IR disc emis- 
sion originates primarily from dust, this conclusion need not 
apply to the gas in the disc. A plausible model might envis- 
age rapid dust agglomeration in the inner disc at the end 
of the CTTS phase, producing the observational effects dis- 
cussed above, while the gas is dispersed on a more leisurely 
timescale. If this were the case, we would expect there to 
be systems where signatures of accretion were visible (e.g. 
strong Ha emission), but which lacked disc emission in the 
infrared. 

If magnetospheric accretion is commonplace, the mod- 
els predict that the colours of T Tauri stars evolve into 
a region of the colour-colour plane characterised by near- 
photospheric K-L colour, and disc-like K-N colour. This oc- 
curs for all values of the stellar magnetic field strong enough 
to disrupt the inner disc at CTTS accretion rates. Strikingly, 
no systems are observed in this region of the diagram. Nu- 
merically, our model tracks move into this 'empty' region 
of the diagram at accretion rates of M ~ 10 -8 Mr^yr" 1 . 
This suggests either that discs are dissipated prior to stars 
reaching this stage (which would imply that passive repro- 
cessing discs are likely to be rare), or that our models are 
missing significant extra sources of infrared flux that would 
shift the tracks redward into consistency with the data. The 
latter possibility cannot be ruled out, though if the extra 
flux arises from a smaller magnetosphere (for example be- 
cause the magnetic coupling is weaker than assumed here), 
the implication would be that stellar angular momentum 
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Figure 5. K-L and K-N colours (from Kenyon & Hartmann, 
1995) for the CTTS with measured photometric periods P* in 
the sample of Bouvier et al. (1995). The upper shded area shows 
the theoretical prediction for a magnetic disc model truncated at 
fl m = R c , with M = 3x 10 — 7 Mq yr _1 , viewed at angles between 
face-on and inclined by 60° (other assumptions as stated in the 
text). The lower shaded area shows the predictions for a pure 
reprocessing disc around the same star with M = 0. 



regulation via magnetic coupling to the disc would be very 
difficult to accomplish. 

Since we have already emphasized that there is consid- 
erable uncertainty in our colour calculation, we note that 
the same conclusion (that reprocessing dominated discs ap- 
pear to be rare) follows from noting that there is no apparent 
pile up of systems at the lower boundary of the CTTS colour 
distribution, as would be expected if there were a long-lived 
population of discs with very low accretion rates. 



4 ROTATION-COLOUR CORRELATION 

In the preceding Section we considered the infra-red colour 
distribution of a large sample of T Tauri stars. For a limited 
subset of these, important additional information in the form 
of photometric periods is available. If this also represents the 
stellar spin period it provides a direct measure of the loca- 
tion of corotation, which in some models is in turn approxi- 
mately equal to the magnetospheric radius. More generally, 
if the magnetic coupling between the star and its disc is 
to act to regulate the stellar angular momentum, then the 
rapid dipolar fall-off of B z (R) implies that i? m must lie close 
to R c . This scenario then predicts that R m should increase 
with increasing P„, and lead to a correlation between the 
stellar spin period and the infra-red colour. 

Figure § shows the K-L and K-N colours of the CTTS 
with measured photometric periods in the sample compiled 
by Bouvier et al. (1995). We make the standard assumption 
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that the photometric period traces the rotation rate of the 
stellar surface and plot the colours as a function of P t . The 
sample comprises 22 systems, all of which have measured K- 
L colours in the compilation of Kenyon & Hartmann (1995), 
and all bar one of which have measured K-N colours. 

We make no pretence that this sample is in any way 
complete or unbiased. However the systems sampled here 
do cover a broad span of stellar periods, from 2 to 12 days, 
which represents a factor of about 3 in the implied radius 
of corotation. No trend in the infra-red colours is seen with 
varying stellar rotation period. This is obvious by eye, and a 
K-S test (for example) readily confirms that the K-L colours 
for the sample split into 'slow' and 'fast' rotators are con- 
sistent with these being drawn from a single distribution. 

The Figure also shows the predictions of a steady-state 
version of the theoretical disc model described in Section 2 
surrounding stars with varying rotation rates. We assume 
7? m = J?c, and take for the stellar parameters M* = Mq, 
R, = 2 R e , B, = lkG, and T, = 4400 K. The magnetic 
heating term is included. The hashed area shows the range 
of colours expected for a disc model with an accretion rate 
of 3 x 10 -7 Moyr -1 , viewed at different inclinations to the 
line of sight. Colours for a reprocessing disc devoid of any 
internal source of luminosity are also plotted. 

Experimentation with the models shows that judicious 
variation of M and B* can roughly reproduce the range of 
colours of these CTTS (though the reddest stars with long 
rotation periods would demand very high accretion rates). 
However, if P* is assumed to correlate well with i? m , the 
models all predict a strong dependence of K-L colour with 
rotation period. This correlation is not seen in the current 
data. Note that K-N colour is not expected to vary greatly 
with P, , as most of the 10 /im flux comes from greater radii 
and is unaffected by magnetospheres of this size. 

The observations display a prominent scatter that could 
be due to variations in accretion rate at a given rotation 
period or to factors not modelled in this analysis. However 
the dispersion, although large, would not mask the expected 
trend if it were uncorrelated with P„ . In particular, if we 
assume that the model with M = 3 x 10~ 7 Moyr -1 defines 
the predicted theoretical trend, then the distribution of (K- 
L)dat a - (K-L) mo< jei does differ significantly between slow and 
fast rotators. Better samples and modelling of the colours 
will be required to determine if this remains a problem for 
magnetospheric accretion models. 



5 DISC EVOLUTION IN BINARY SYSTEMS 

A large fraction of T Tauri stars are members of binary or 
multiple systems (Chez 1996), and some of these are close 
enough to influence the evolution of circumstellar discs. We 
consider here the evolution of discs in these moderately close 
(~ 10 — 100 a.u.) binaries. 



5.1 Truncation radii 

Previous theoretical work has shown that the gravitational 
influence of a binary system on circumstellar material is ex- 
pected to be twofold: the individual circumstellar discs will 
be truncated by the tidal effects of the companion star and 
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Figure 6. The assumed outer radius of the accretion disc of the 
primary (upper curve) and secondary (lower curve) as a fraction 
of the binary separation, for various mass ratios q. Adapted from 
Papaloizou & Pringle (1977). 



a gap will be cleared between the binary and any surround- 
ing material forming a circumbinary disc. The truncation 
and gap-clearing processes occur rapidly, on a dynamical 
timescale, and so for our purposes the initial conditions for 
discs in a binary can be taken as two truncated circumstel- 
lar discs. These discs will then evolve independently, un- 
less a circumbinary disc exists and material is able to flow 
across the gap (for a review, see Lubow & Artymowicz 1996). 
This process is predicted to occur, though observationally 
circumbinary discs do not appear to be common (Mathieu 
1996). 



5.2 Discs in systems of varying separation 

The simplest effect of binary companions on disc evolu- 
tion might be expected to be a systematic variation of disc 
lifetime with binary separation. The reduced viscous time 
tu — -Rout/^out at the outer edge should lead to faster disc 
evolution in closer systems, and a higher fraction of WTTS 
as compared to CTTS. This effect has been searched for 
(Chez, Neugebauer & Matthew 1993; Simon & Prato 1995), 
but although there are systems such as DI Tauri (Meyer 
et al. 1997) that seem to show the expected effect there 
appears to be no fully convincing evidence for a general 
trend of CTTS fraction with separation. Extension of the 
observations to smaller separation binaries should be able 
to determine if this implies ubiquitous replenishment of cir- 
cumstellar discs from a circumbinary reservoir, as is already 
implied from the mere existence of discs in the very closest 
(spectroscopic) binaries. 



5.3 Discs in systems with unequal mass ratios 

For binaries with unequal mass ratios, we can also ask 
whether the different truncation radii of discs surrounding 
the primary and secondary will lead to observable differen- 
tial disc evolution. Figure shows the predicted outer radii 
of tidally truncated discs in a binary system as a function 
of the mass ratio q = M secon dary /Mprimary (results taken 
from Papaloizou & Pringle 1977; see also Paczyhski 1977). 
We take the 'primary' to be the more massive star, so that 
q < 1. For binaries with q ~ 0.5, the ratio of the outer radii 
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Table 1. Summary of the initial conditions for the binary calcu- 
lations 
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of their accretion discs is expected to be ~ 1.5, rising to 
2 — 3 for more extreme systems with q ~ 0.1. This dispar- 
ity in radii will affect the evolutionary timescales of the two 
discs, while the epoch at which the transition from CTTS 
to WTTS status occurs will also depend on the disc initial 
conditions. 

If the discs survive for long enough the accelerated evo- 
lution of the secondary's smaller disc compared to the pri- 
mary's will eventually overwhelm any differences in the ini- 
tial conditions of the two discs. Whatever the initial ratios of 
surface density, disc mass and accretion rate, at late enough 
epochs we thus expect that the secondary should display the 
weaker disc of the pair. Our numerical results aim to quan- 
tify when that switchover should occur, and how the models 
compare to the growing body of observations of close binary 
systems. 

Table [j] summarises the main parameters for the numer- 
ical runs. In all cases we consider a binary with separation 
50 a.u., take a = 5x 10 -3 , and impose a, vr — boundary 
condition at Rout- The smaller value of a gives similar disc 
mass fractions as for the models in Section 3. 600 radial grid 
points are used. 

The initial conditions for the primary are taken as 
a steady-state disc with an accretion rate of M = 3 x 
10 -7 MqJT -1 , that fills its Roche lobe out to R ou t- For the 
q — 0.44 case a weakened initial disc was also run, with the 
same initial accretion rate but with the radial extent of the 
disc truncated such that Mdi sc was reduced by a factor of 
approximately four. 

For the discs surrounding the secondaries, our stan- 
dard case commences with an accretion rate equal to that 
of the primary. These discs are initially a factor 2-6 times 
less massive than the corresponding primary disc. However, 
cognisant of the fact that in some formation scenarios the 
secondary disc might accrete more infalling matter than the 
primary (Bate & Bonnell 1997), we also consider initially 
higher accretion rates through the secondary's disc. These 
are chosen such that the Q parameter, 

Q = -7^ (10) 

equals unity at R ou t, and represent the maximum plausible 
accretion rates through steady discs described by the verti- 
cally averaged equations. For Q < 1 gravitational instability 
sets in (Laughlin & Bodenheimer 1994), and leads to a rapid 




Figure 7. The decline of the disc mass, surface density at 1 a.u., 
and accretion rate (in M0yr — 1 ), for the q = 0.44 calculation. The 
solid curves show the primary disc for the standard (upper curves) 
and low mass (lower curve) disc initial conditions; the dashed lines 
show the secondary disc for high and low initial accretion rates. 
Here we have assumed a binary separation of 50 a.u., and taken 
a = 5 X 10" 3 . 



redistribution of disc angular momentum. A more massive 
disc would thus be expected to relax to a gravitationally 
stable configuration on a timescale that is short compared 
to the viscous timescales over which we compute the evolu- 
tion. These high M secondary discs have surface densities 
that are initially considerably higher than for the primary, 
and masses that are comparable, or (in the q — 0.44 case) 
slightly larger. 

Figure ^ shows in detail the evolution of the discs for 
the runs intended to mimic a q = 0.44 binary system. Three 
possible indicators of the strength of the circumstellar discs 
are plotted, the total disc mass, the surface density at a 
radius of 1 a.u., and the accretion rate, which is expected 
to correlate well with measures of activity based on Ha and 
UV flux. 

For the standard initial conditions of the primary disc, 
and for all three measures of disc strength, the faster evolu- 
tion of the smaller disc around the secondary is sufficient to 
weaken that disc below that of the primary in a relatively 
short time ~ 10 6 yr or less. The initial conditions matter 
most for measures based on the accretion rate, but even 
here the primary disc is substantially stronger than that of 
the secondary at accretion rates well above the typical de- 
tectable threshold. 

If the disc of the primary is initially substantially 
smaller than the tidal truncation radius then a much longer 
period is required before the secondary's disc becomes 
weaker than that of the primary. In particular, the accre- 
tion rate through the primary disc remains below that of the 
models for the secondary down to M ~ 3 x 10~ 9 M®yr , 
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which would be hard to measure unambiguously observa- 
tionally. Hence if the initial conditions were as biased in 
favour of a strong secondary disc as in this calculation they 
would overwhelm the faster evolution of the secondary's disc 
for all observations dependent on the disc accretion rate. 

5.4 Variation with mass ratio 

Figure |§] summarises the results of the runs in Table |l|. For 
each mass ratio and indicator of disc strength (disc mass, 
accretion rate, and surface density at an arbitrary radius 
of 1 a.u.) we plot the relative strength of the primary and 
secondary discs as a function of time. As in Section 3 we 
caution that the absolute time scaling of the models is some- 
what arbitrary - a has been chosen so as to give plausible 
disc lifetimes in the Myr range. 

As we have already discussed, for q = 0.44 the evo- 
lutionary timescales of the two discs are sufficiently simi- 
lar that the initial conditions are important through to late 
times. For reasonable choices, there is a lengthy period when 
the primary and secondary discs have at least comparable 
strength. Observations of such binaries thus potentially re- 
tain information on the formation histories of the discs. 

Few binaries appear to have very extreme mass ratios. 
However studies of those systems can provide a stringent test 
of the evolution of the disc, since for q < 0.2 — 0.3 the large 
disparity in evolutionary timescales is sufficient to rapidly 
overwhelm variations in the initial disc masses and radii. 
From Fig. |^ it can be seen that for the two most extreme 
mass ratios essentially all binaries are predicted to have a 
stronger circumprimary disc than circumsecondary. This ap- 
plies for all measures of disc strength, though the disc mass 
is always the most robust measure. 

The results can also be used to estimate the predicted 
fraction of 'mixed' systems, i.e. the number of binaries pair- 
ing a WTTS with a CTTS as a fraction of the total num- 
ber of binaries containing at least one CTTS. To do this 
we assume that the transition from CTTS to WTTS occurs 
at some fixed value of the disc mass (or accretion rate or 
surface density) for both primary and secondary, and use 
the numerical models to provide the time spent as a doubly 
strong (CTTS + CTTS) and as a mixed (CTTS + WTTS) 
system. 

For the q — 0.44 run, the results depend strongly on 
when the epoch of disc dispersal occurs. If the disc is lost rel- 
atively early on, when only a few viscous times have elapsed, 
then the predicted fraction of mixed systems depends almost 
entirely on the uncertain initial conditions. 

For lower-mass-ratio systems, where the disparity in 
evolutionary timescales is greater, the influence of the initial 
conditions is less important. For q — 0.19 the time spent in 
the 'mixed' state is expected to exceed that spent as a dou- 
bly strong system, while for q = 0.1 the ratio of times is 
~ 5 for all measures of activity (i.e. the fraction of mixed 
systems amongst all binaries containing at least one CTTS 
would be ~ 0.8). This applies even if the discs are destroyed 
at an accretion rate of order 10 -8 M0yr _1 . The observation 
of even a few systems where the mass ratio was this extreme 
can thus provide a stringent test - the standard models pre- 
sented here are incompatible both with a low fraction of 
mixed (CTTS, WTTS) systems, and with systems where 
the primary is weaker than the secondary at any except the 



highest accretion rates. Observations to the contrary would 
be good evidence for accretion from a circumbinary disc. 

5.5 Discussion 

In this Section we have investigated the evolution of discs 
in close binary systems. We find that the secondary's disc 
should always become weaker than that surrounding the pri- 
mary at a sufficiently late epoch, but that the timing of the 
switchover is a strong function of the mass ratio of the bi- 
nary. For mass ratios of ii 0.5, a secondary that has an 
initially higher accretion rate through its disc might remain 
stronger than the primary until the disc is cleared. For more 
extreme mass ratios, we expect the faster evolution of the 
smaller secondary disc to dominate over the influence of the 
initial conditions. At late times the secondary should display 
weaker disc signatures than the primary, and there should be 
a relatively high fraction of mixed T Tauri binaries pairing 
a CTTS primary with a WTTS secondary. 

Current observations appear to be broadly consistent 
with this picture. Brandner & Zinnecker (1997) studied a 
sample of binaries in the 90-250 a.u. range, finding 3 sys- 
tems (out of 12 containing a CTTS) where a CTTS sec- 
ondary was paired with a WTTS primary. In this sample, 
all the measured mass ratios were q > 0.5, and the analy- 
sis was based on Ha flux. This sample thus falls into the 
category where we would predict that the initial conditions 
of the discs formed around the stars should be important, 
and strong accretion activity among the secondaries is not 
unexpected. 

The photometric properties of close T Tauri binaries 
have been studied by Prato & Simon (1997), and Chez, 
White & Simon (1997). These studies show that mixed 
(CTTS, WTTS) binaries appear to be rare - none were 
found in the sample of 12 systems examined by Prato & 
Simon. This is probably marginally consistent with our mod- 
els, at least if discs in these binary systems are cleared at a 
relatively high accretion rate. However, it is clear that a sim- 
ilar observational result for systems with an extreme mass 
ratio would be inconsistent with these independent evolution 
models, and would point either to common replenishment of 
the discs from a circumbinary disc, a much weaker depen- 
dence of t v on R than that assumed here, or a coordinated 
mechanism for their common destruction. 



6 CONCLUSIONS 

In this paper we have presented models for the evolution 
of circumstellar discs around T Tauri stars. The models 
combine a pre-main-sequence stellar evolution track with 
a model for viscous disc evolution. Over the estimated 1-10 
Myr lifetime of T Tauri discs (Strom 1995) there are large 
changes in both the stellar luminosity and the disc accretion 
rate (Hartmann et al. 1998), making the inclusion of both 
components essential in an evolutionary model. 

Comparison of our theoretical tracks with the distribu- 
tion of infra-red colours of T Tauri stars in Taurus- Auriga 
(Kenyon & Hartmann 1995) yields two main conclusions. 
Most securely, in all models the transition in the colour- 
colour plane between systems looking like Classical T Tauri 
stars, and those appearing as Weak-lined systems, is slow, 
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Figure 8. The relative strength of the primary and secondary discs, as a function of time, for varying mass ratios q. Disc 'strength' is 
measured by disc mass (left panels), accretion rate onto the star (centre panels), and surface density at a radius of 1 a.u. (right panels). 
The solid lines show the ratio of the strength of the primary to secondary discs for models where the secondary disc initially has a high 
accretion rate. The dashed lines show the case where the two stars have initially identical accretion rates. For a mass ratio q = 0.44 we 
also show with the dotted line the case where the primary disc is initially of low mass, as described in the text. 



and occurs on the viscous timescale of the outer disc. This 
occurs even if a dynamically significant stellar magnetic field 
disrupts the inner disc. However observationally this tran- 
sition must occur very rapidly, since there are few systems 
with properties intermediate between CTTS and WTTS. 
Non-viscous processes must be the agents of this transition. 

A more detailed comparison of the models to the data 
provides some constraints on when disc dispersal must oc- 
cur. The presence of dynamically significant stellar magnetic 
leads to a distinctive evolutionary track of stars in the (K-L)- 
(K-N) plane. As the accretion rate declines, the (K-L) colour 



drops rapidly to close to bare photospheric levels while the 
system is still red in (K-N). No systems with this combina- 
tion of colours are observed, which implies that the models 
can only be consistent with the observations if purely repro- 
cessing discs are rare. We speculate that changes in the disc 
structure at the epoch when the disc goes from being ac- 
tive (heating dominated by viscous processes), to passively 
reprocessing stellar radiation, could be responsible for initi- 
ating the disc's rapid demise. 

Applying the models to close binary systems, we find 
that the observed lack of systems pairing CTTS with WTTS 
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is consistent with disc destruction occurring at a similar ac- 
cretion rate to that inferred from the single star models. The 
models suggest that the frequency of mixed T Tauri binaries 
should increase greatly at extreme mass ratios. Observations 
of such systems, and studies that probe the lifetime of the 
CTTS phase as a function of separation, are a good indica- 
tor of T Tauri disc evolution and the ubiquity of accretion 
from circumbinary discs. Conversely, observations of discs 
in binaries with more nearly equal mass components probe 
the formation history of those discs, especially if observed 
at early epochs. 

For the subsample of systems with measured photo- 
metric periods, we find no correlation between the infra- 
red colours and the photometric period, which we assume 
is probably a good tracer of the stellar rotation rate. If the 
magnetospheric radius was approximately equal to the coro- 
tation radius, or a fixed large fraction of it, a correlation of 
(K-L) with P* would be expected. The lack of such a cor- 
relation, along with suggestions that the magnetosphere lies 
well inside the radius of corotation (Kenyon, Yi & Hartmann 
1996; Meyer, Calvet & Hillenbrand 1997), pose potentially 
serious problems for models that seek to explain the slow 
rotation of Classical T Tauri stars as a result of magnetic 
linkage between the star and its disc. 
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